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Abstract 

A capability to calculate surface heating rates has 
been incorporated in an approximate three-dimensional 
inviscid technique. Surface streamlines are calculated 
from the inviscid solutionr and the axisymmetric analog 
is then used along with a set of approximate convective- 
heating equations to compute the surface heat trans- 
fer. The method is applied to blunted axisymmetric and 
three-dimensional ellipsoidal cones at angle of attack for 
the laminar flow of a perfect gas. The method is also ap- 
plicable to turbulent and equilibrium-air conditions. The 
present technique predicts surface heating rates that com- 
pare favorably with experimental (ground-test and flight) 
data and numerical solutions of the Navier-Stokes (NS) 
and viscous shock-layer (VSL) equations. The new tech- 
nique represents a significant improvement over current 
engineering aerothermal methods with only a modest in- 
crease in computational effort. 


Nomenclature 


A, B, V, J 

geometric factors 


tangential unit vectors on body surface 

Gj; 7 6 r , 

unit vectors of cylindrical coordinate system 


unit vectors of shock-oriented coordinate 
system 

e £ ; e /3 ; e n 

unit vectors of streamline coordinate system 

/ 

shock radius 

/ 

body radius 

h(,hp 

scale factors of shock-oriented coordinate 
system 

h(, hp 

scale factors of streamline coordinate system 
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M Mach number 

n coordinate normal to shock 

n coordinate normal to body 

p static pressure 

q heat-transfer rate 

R radius of curvature 

и, v,w velocity components of shock-oriented 

coordinate system 

V velocity magnitude 

V velocity vector 

x,r,<f) cylindrical coordinate system 

x,y,z Cartesian coordinate system 

a angle of attack 

T shock angle relative to freestream velocity 

T body angle relative to freestream velocity 

8 p shock angle in circumferential direction 

8 $ body angle in circumferential direction 

rj stream function ratior’f' /\? s 

9 inclination angle of surface streamlines 

к, £,K /3 shock curvatures 

£,/? shock coordinates 

£,/? streamline coordinates 

p density 

a shock angler^ — 8$ 

<7 body angler^ — 8 $ 

stream functions 

Subscripts 

b body 

s shock 

w wall 

oo freestream conditions 


Introduction 

The thermal design of hypersonic vehicles involves 
accurately and reliably predicting the convective heat- 
ing over the surface of the vehicle. Such results may 
be obtained by numerically solving the Navier-Stokes 
(NS) equations! 1 or one of their various subsets such as 
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the parabolized Navier-Stokes (PNS) 2 and viscous shock- 
layer (VSL) equations 3,4 for the flowfield surrounding the 
vehicle. HoweveiT due to the excessive computer stor- 
age requirements and run times of these detailed ap- 
proachesrthey are impractical for the preliminary design 
environment where a range of geometries and flow pa- 
rameters are to be studied. On the other handrengineer- 
ing inviscid- viscous methods 5-8 have been demonstrated 
to adequately predict the heating over a wide range of 
geometries and aerothermal environments. Various ap- 
proximations in the inviscid and boundary-layer regions 
reduce the computer time needed to generate a solution. 
This reduction in computer time makes the engineering 
aerothermal methods ideal for parametric studies. 

Two of the simpler engineering aerodynamic heating 
methods that are currently used are AEROHEAT 5,0 and 
INCHES.' Both use the axisymmetric analog concept 9 
which allows axisymmetric boundary-layer techniques to 
be applied to three-dimensional (3-D) flows provided the 
surface streamlines are known. AEROHEAT calculates 
approximate surface streamlines based solely on the body 
geometry. INCHES uses an approximate expression for 
the scale factor in the windward and leeward planes which 
describes the spreading of surface streamlines. Circum- 
ferential heating rates are then generated by an empirical 
relation. Another area of approximation is the surface 
pressure distribution employed by the engineering meth- 
ods. AEROHEAT assumes modified Newtonian theory 
which is inaccurate for slender bodiesr while INCHES 
uses an axisymmetric Maslen technique. 10 The deficien- 
cies and limitations of these approximations to the surface 
streamlines and pressures in the engineering aerothermal 
met.hodsr along with their corresponding effects on the 
surface heat transfeiThave been documented in Refs. 11 
to 13. 

An approximate 3-D inviscid method 14,15 has been 
developed that is more accurate than modified Newto- 
nian theory and has a wider range of applicability than 
the axisymmetric Maslen technique. The inviscid tech- 
nique uses two stream functions that approximate the 
actual stream surfaces in tint shock layer and a modified 
form of the Maslen second-order pressure equation. 10 The 
method has been shown to calculate the inviscid flowfield 
about. 3-D blunted noses as well as 3-D afterbodies reason- 
ably accurately and much faster than numerical solutions 
of the inviscid (Euler) equations. 14 

In this papeiTt.he approximate inviscid technique em- 
ploys the axisymmetric analog to predict laminar and 
turbulent surface heating rates using the approximate 
convective-heating equations of Zoby et. al . 1 ' Both perfect 
gas and equilibrium-air flows are considered. Improved 
>n r face streamlines are calculated based on both the body 
geometry and surface pressure distribution. Surface heat- 
ing rates are presented for spherically-blunted and asym- 
metric ellipsoidal cones at angle of attack. Comparisons 
are made between results of the present techniquerVSL 
and NS solut.ionsr and available experimental data to 



Figure 1. Shock wave geometry: side view. 

demonstrate the accuracy and capability of the present 
engineering technique. 

Analysis 

This section describes the 3-D inviscid techniquerthe 
procedure for computing inviscid surface streamlinesEand 
the application of the axisymmetric analog. Approxima- 
tions and coupling issues are also discussed. 

Inviscid Method 

Since a detailed description of the approximate 3-D 
inviscid method has been presented previouslyl 14,15 only 
a brief outline of the inviscid method is given here. 

Coordinate Systems 

The three-dimensional shock surface can be repre- 
sented by 

r s = f(x,ft) (1) 

where (x,r,<f>) are wind-oriented cylindrical coordinates 
with corresponding unit vectors (e F , e r , e^,). The r-axis is 
aligned with the freestream velocity vector and is normal 
to the shock surface at the origin. Two anglesE^r, ft) 
and r(r, c/))r describe the shock wave shape and are de- 
fined as 

, c 1 d f , p c m 

tan 8,p = -— tan T = — cos ( 2 ) 

An additional angle is given by a = ft — 6$. All angles 
are shown in Figs. 1 and 2. For the special case of ax- 
isymmetric flowr?’ s = f(x)TT = r(r)rii^ = OTarid a = ft. 

NextTa. shock-oriented curvilinear coordinate system 
(^, ft, n) is defined where ^ and ft represent coordinates of 
a point on the shock surface and n is the inward distance 
normal to the shock. Differential arc lengths along each 
coordinate direction at the shock are hp dftT hp d/3T and 
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at the shock. Along a line normal to the shockr these 
expressions are 



Figure 2. Shock wave geometry: rear view. 


dn where h p and hp are scale factors for the correspond- 
ing coordinates. This coordinate system is well-suited for 
hypersonic flow (M m 1) and thin shock layers. 

The unit vectorsFe^ and e^Tare tangent to the shock 
surface and are chosen such that is in the direction of 
the tangential velocity just inside the -shock surface. The 
unit vector ep is then defined to be perpendicular to 
and e n . In cylindrical coordinatesrthe unit vectors of the 
curvilinear coordinate system are given by 

= cosr e F + sinr (cosi 5p e F — siin^ e^,) 
ep = sindi^ e r + cosd^ ep (3) 

e n = sinr e F — cosT (cosi 5p e F — sirnS^, epft 

Although this curvilinear coordinate system is orthogo- 
nal at the shock surfacer it is nonorthogonal within, the 
shock layer for a general three-dimensional shock. How- 
eveiTfor thin shock layersrorthogonality may be assumed 
everywhere. 

The velocity is defined in terms of the unit vectors 
at the shock as 


V = we^ + ve n + wep (4) 

From the definition of and e^Tthe crossflow velocity 
component at the shockr®,, Ids equal to zero. 

Governing Equations 

The governing equations for 3-D inviscid flow are 
simplified by assuming that the velocity component w 
is equal to zero not only at the shock but throughout the 
shock layer. This yields two stream functionsr <F (which 
is equal to ft here) and 'FT which approximate the actual 
stream surfaces in the shock layer. The stream function 
’F is analogous to the Stokes stream function for axisyrn- 
rnetric flow. 

Approximate expressions for the pressure and nor- 
mal velocity component are then obtained by transform- 
ing the normal momentum and continuity equations to 
streamline coordinates and evaluating the flow variables 


p(V ) = Ps +Pi (V ~ 1) +P2 (if ~ 1) (5) 

v(r)) = v s + vi (r) - 1) (6) 


where 


and 


Pi = 
P2 = 
v l = 


hp 

\F s v s tanr 

2 hp 

hp cosr 


( K e + k p ) 


(Kp+Kp) 


h = 


^7 


Defining 4t = 0 to be the body surface gives i] = 1 on 
the shock and i] = 0 on the body. Note that Eq. (5) 
reduces to Maslen’s second-order pressure equation 113 for 
axisymmetric flow if the scale factor hp is equal to the 
shock radius r s . 


Method of Solution 

Since the inviscid method is an inverse onerthe shock 
shape must be varied until the correct body shape is pro- 
duced. The resulting iteration procedure is handled dif- 
ferently in each region of the flow. 

In the stagnation region of a blunt body traveling 
at hypersonic speedsIThe flow is subsonic and the shock 
shape for the entire subsonic-transonic region must be 
determined globally. A 3-D shock given by longitudi- 
nal conic sections blended in the circumferential direction 
with an ellipse is assumed. The parameters describing the 
shock are iterated until the body shape (4t = 0) generated 
by the approximate inviscid method matches the actual 
body shape at several discrete points. In this st.udyFsix 
shock parameters are varied until the calculated body is 
matched to the actual body at six locations. 

Once past the transonic regionr the inviscid flow is 
totally supersonic and a marching scheme is well posed. 
The shock surface from the transonic region forms a start- 
ing solution for the marching procedure. The shock vari- 
ables are extrapolated in £ along a number of constant 
ft lines which circle the shock. On each liner the shock 
curvature up is locally iterated until the calculated body 
shape matches the correct body. The shock variables are 
then advanced downstream to the next ^-location and the 
process repeated. 


Axisymmetric Analog 

The 3-D boundary-layer analysis is simplified by us- 
ing the axisymmetric analog 9 as is done in most engi- 
neering aerothermal methods. The 3-D boundary-layer 
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equations are first written in a streamline coordinate sys- 
tem. The crossflow velocity component tangent to the 
surface but normal to the streamline is then assumed to 
be zero. This simplification reduces the 3-D boundary- 
layer equations to the axisymmetric form provided the 
distance along the streamline is substituted for the surface 
distance and the scale factor describing the divergence of 
the streamlines is interpreted as the axisymmetric body 
radius. Axisymmetric boundary-layer methods can then 
be employed in the existing 3-D inviscid technique. 


The orientation of the inviscid surface streamlinesr 
given by 6Tis found by applying the momentum equations 
along the body surface using the pressure distribution 
generated by the inviscid solution. By writing the mo- 
mentum equations in streamline coordinatesr taking the 
scalar product with e^Tand substituting the unit vectorsr 
Eqs. (9) and (lO)Tthis may be expressed as 

1 89 sinf da 11 dpb 
hp dp h p dp PbVp hp dp 


Inviscid Surface Streamlines 

Before applying the axisymmetric analogr inviscid 
surface streamlines are computed from the approximate 
inviscid solution. Inviscid surface streamlines may be 
calculated from the surface pressure distribution 5 or 
from the velocity components. 8 The approximate inviscid 
method 14,15 used here predicts accurate surface pressures! 1 
but the direction of the velocity on the surface is not ac- 
curate. Thereforerin the present methodrstreamlines are 
calculated from the surface pressures. 

A streamline coordinate system 5 (p,P, n) is defined 
where f and [3 are coordinates of a point on the body 
surface and h is the distance normal to the body. The 
bars indicate the variables apply to the body and not 
the shock. Differential arc lengths along each coordinate 
direction at the body are hp d^Thp dpT and dh where h p 
and hp are scale factors for the corresponding coordinates. 
If the body surface is represented by rj = f(x, p) in wind 
axes with the axial coordinate parallel to the freestream 
velocity and passing through the stagnation pointT the 
unit vector normal (outward) to the body surface is given 

by 


e„ = — sinr e x + cosT (cos8p e r — sini)^ ep) (7) 

The body angles are defined in the same fashion as the 
shock angles and are 

, 7 1 df f 8f - 

taniu = - tan I = — cos bp (8) 

/ 8(p 8x 

The tangential unit vectors at the surfaceFeg and e^Tare 
similar to the tangential unit vectors at the shock. From 
Ref. 5Tthey are given as 

e p = cos 9 e y + sin# ej (9) 

ep = — sin# eg + cos# ef (10) 

where 

eg = cosT e x + sinr (cosbp e r — sini)^ ep) (11) 

ef = sini)^ e r + cosi 5p ep (12) 

and the angle 9 represents the orientation of the surface 
streamlines. Note that the vectorsFeg and efTare identi- 
cal in form to the unit vectorsFeg and e^Tdehned at the 
shock. 


where d = if) — bp. The scale factor hp can be determined 
by noting that for an orthogonal curvilinear coordinate 
system 

m <*»•»> = w 

Taking the scalar product of this equation with and 
again substituting the unit vectorsr Eqs. (9) and (10)T 
yields 

1 8 In ha 1 89 sinf 8d 

=-£- = -| (14) 

% hp dp hp dp 

Equations (13) and (14) may be integrated along a 
surface streamline to obtain the streamline direction 9 
and the scale factor hp. Although the surface stream- 
lines can be determined after the inviscid solution has 
already been calculatedr it was found to be more con- 
venient to compute the inviscid solution and the surface 
streamlines simultaneously. Before applying these equa- 
tions along shock coordinatesr transformation operators 
relating derivatives with respect to the the streamline co- 
ordinates (£,,P) to derivatives with respect to the shock 
coordinates (£,/?) are needed. In the approximate invis- 
cid methodrthe curvilinear coordinate system is assumed 
to be orthogonal throughout the shock layer. This as- 
sumption simplifies the analysis but does not change the 
form of the approximate pressure and velocity relationsr 
Eqs. (5) and (6)rsince the flowHeld variables are evaluated 
at the shock where the coordinate system is orthogonal. 
HoweverFat the body surfacerthe correct coordinate di- 
rections need to be considered. Following the approach 
of Ref. 15 and using the nonorthogonal directions at the 
surfacerthe transformation operators are 

= (Be r e ( -Ve r e p )^--^ 

+ (-1 >e f - e f+ .4e r e,)i-fl (15) 

= (Bep-e^-Vep-ep)^--^ 

+ (~Vep -e^ + Aep -ep) (16) 


and 


3_d_ 

hpd^ 


3_d_ 

hp dp 


where 


A = 1 — njKf 
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10.6 Present method 

1.1 in Present method O Reentry F (Ref. 22) 

O Exp. (Ref. 21) 
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